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Context. Photospheric absorption lines in classical T Tauri stars (CTTS) are weak compared to normal stars. This so-called veiling is 
normally identified with an excess continuous emission formed in shock-heated gas at the stellar surface below the accretion streams. 
Aims. We have selected four stars (RW Aur A, RU Lup, S CrA NW and S CrA SE) with unusually strong veiling to make a detailed 
investigation of veiling versus stellar brightness and emission line strengths for comparisons to standard accretion models. 
Methods. We have monitored the stars photometrically and spectroscopically at several epochs . 

Results. In standard accretion models a variable accretion rate will lead to a variable excess emission. Consequently, the stellar 
brightness should vary accordingly. We find that the veiling of absorption lines in these stars is strongly variable and usually so 
large that it would require the release of several stellar luminosities of potential energy. At states of very large line dilution, the 
correspondingly large veiling factors derived correlate only weakly with brightness. Moreover, the emission line strengths violate the 
expected trend of veiling versus line strength. The veiling can change dramatically in one night, and is not correlated with the phase 
of the rotation periods found for two stars. 

Conclusions. We show that in at least three of the stars, when the veiling becomes high, the photospheric lines become filled-in by 
line emission, which produces large veiling factors unrelated to changes in any continuous emission from shocked regions. We also 
consider to what extent extinction by dust and electron scattering in the accretion stream may affect veiling measures in CTTS. We 
conclude that the degree of veiling cannot be used as a measure of accretion rates in CTTS with rich emission line spectra. 
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1. Introduction 

Classical T Tauri stars (CTTS) are low-mass pre-main-sequence 
objects with intense emission line spectra. In standard models 
the spectacular emission spectrum is explained as the manifes- 
tation of magnetospheric accretion of material from a circum- 
stellar disk. Gas is heated while accelerated towards the stel- 
lar surface, where energy dissipates in strong shocks. Since 
the stars may host large-scale organized magnetic fields the 
shocked gas, producing an excess continuous emission, would 
be distributed over areas located at high latitudes surrounding 
the magnetic poles. In addition, the stars show flares and chro- 
mospheric/coronal/wind activity. For reviews of properties and 
models of CTTS see Petrov d20Tj3l i and Bouvier et al. {2007 a) . 

The photospheric absorption lines are diluted (an effect 
known as veiling). By comparing such spectra with those of non- 
veiled stars one can derive the energy distribution of the veiling. 
These distributions have been matched with Paschen continuous 
emission as well as black body radiation, and for CTTS the de- 
rived temperatures are of the order of 10 4 K (see e.g. Basri & 
Batalha [1^901 Hart igan e t al. [1991] [19951 Val enti et al. [19931 
Calvet & Gullbring [T998l Gullbring et al. [19981 . The amount of 
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veiling has been used to estimate accretion rates in TTS. One 
expects that when the accretion rate changes, the veiling and the 
stellar brightness also will change accordingly. However, little 
systematic, long-term spectroscopic and photometric monitor- 
ing of TTS has been performed to explore this prediction. The 
first studies of this kind, based on low- to medium-resolution 
spectroscopy, gave the confusing result that although some stars 
appeared to obey the expected correlation between the degree of 
veiling and stellar brightness, other stars, and especially some 
heavily veiled stars, did not (Gahm et al. 119951 Hessman & 
Guenther 1997 ; Chelli et al. 1999 ). Later high-resolution studies 
showed that the predicted veiling-brightness relation does not al- 
ways hold in RW Aur A (Petrov et al. 2001; Petrov & Kozack 
120071 1, while in AA Tau it does to some extent (Bouvier et al. 
2003 2007b). In the present Letter we investigate the behaviour 
of veiling in relation to stellar brightness and emission line fluxes 
in RW Aur A, RU Lup and the binary S CrA (NW and SE). This 
project is part of a monitoring campaign of CTTS involving op- 
tical and NIR spectroscopy, photometry and X-ray observations. 



2. Observations 

High-resolution optical spectra of RU Lup and S CrA were ob- 
tained during consecutive nights in 2000, 2002 and 2005 with 
the UVES spectrograph on the 8-m VLT/UT2 of the European 
Southern Observatory. For descriptions of the observations and 
data reductions see Stempels et al. (2007, Paper 1). The individ- 
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ual spectra of the S CrA components, separated by only 1'.'3, are 
resolved. RW Aur A was observed during several nights each 
year from 1995 to 1999 using the SOFIN spectrograph at the 
Nordic Optical Telescope, and the observations and reductions 
were described by Petrov et al. (2001 , Paper 2). 

We have compiled a long-term record of RU Lup and the two 
S CrA components with UBVRI/JHK photometry since 2003, 
using the ANDICAM dual-channel imager on the SMARTS 1.3 
m telescope at Cerro Tololo. For RW Aur A we derived V mag- 
nitudes from stars in the CCD-field at NOT with exposures taken 
before and after the spectroscopic observations. 



3. Results 

All four program stars are "extreme" CTTS, with very prominent 
emission line spectra and with only weak traces of photospheric 
absorption lines, implying that the veiling is high. We have ob- 
tained precise measures of the veiling factor (VF), defined by the 
ratio of the excess continuum to the intrinsic photospheric con- 
tinuum flux. The VFs were determined by comparing observed 
spectra with synthetic template spectra and, in the case of RW 
Aur A, with observed template spectra. The procedures are de- 
scribed by Stempels & Piskunov (2003), where also the veiling 
of RU Lup is discussed, as well as in Papers 1 and 2. The veil- 
ing was measured in spectral windows containing both visible 
and isolated absorption lines. This allowed us to derive the en- 
ergy distributions of the veiling showing, as is normal, a rise 
towards the blue spectral region. The veiling is strongly variable 
in all stars. No star ever fell below VF = 1 in the V band region 
(corresponding to a 50% decrease in line strength), and values 
exceeding VF =10 were measured. If such high veilings result 
from excess continuous emission, then the flux in V should be- 
come very large, corresponding to several stellar luminosities 
being liberated at the stellar surface from gravitational energy 
deposited by the accretion streams. 

3.1. Veiling versus stellar brightness 

RU Lup shows the same photometric behaviour as in numerous 
earlier investigations (see e.g. Giovannelli 1994 and Pojmanski 
2002). The star rarely varies by more than +0.4'" in V and 
+0. 1 5'" in /. SMARTS data collected nearly simultaneously with 
our spectroscopic observations in 2005 is shown in Fig. Q] (top 
panels) together with veiling factors for the V band. During the 
first two spectroscopic nights the veiling was extremely high (VF 
> 5, peaking at 7.6) after which it declined to more 'normal' 
levels of VF « 2.6. We recorded a small decline in brightness 
of * 0.2™ in all bands over this period. If the changes in veil- 
ing were related to a hot area below the accretion stream, then 
the expected decrease in V should amount to » 1'". Similarly for 
RW Aur A (see Fig.Q] bottom panels), the veiling can decrease 
from VF = 7.0 to 1.5 over 2 days, while the flux in V decreases 
by « 0.1'" compared to 1.3'" expected from a fading hot area. 

Both S CrA components vary dramatically in brightness 
(Walter & Miner 2005), and we recorded large changes in the 
veiling, which on occasion rises to VF > 10. Photometry and 
spectroscopy secured in 2005 show that the coupling between 
veiling and brightness is weak in both components. During this 
run the VF increased steadily from 5.9 to 10.7 in S CrA NW. 
This rise was accompanied by a brightening of only 0.1 m in V, 
while 0.6'" would be expected from a changing hot area. In S 
CrA SE, when VF rose by 30% from Aug. 16 to 17, there was 
no detectable change in VRI, but the star became brighter in B. 
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Fig. 1. Top panels: Overview of veiling factors derived from 
UVES spectra collected over 4 nights of RU Lup in August 2005 
including two nights of extremely high veiling, and SMARTS 
BVRI photometry for the same period. Bottom left: Spectral re- 
gion of RW Aur A obtained at high and low degrees of veil- 
ing. Bottom right: Veiling factor and V magnitude of RW Aur A 
during an event where the veiling declined drastically from an 
extremely high level (circles: our NOT photometry; squares: 
Grankin et al. (I2007T ). 

Hence, in all four program stars, the veiling is unrelated, or 
only weakly related to brightness during states of high veiling. 
These results are inconsistent with the interpretation that veiling 
as measured from high-resolution spectroscopy is related to the 
flux level in the accretion continuum. 

For spectra of RU Lup taken in 2002 we could use a sec- 
ondary flux standard to flux-calibrate our spectra on three nights. 
This allowed us to calculate synthetic B and V magnitudes of 
RU Lup, and directly compare the stellar brightness with veil- 
ing. Individual measures of magnitude and colour are uncertain 
but during one night, when VF increased from 2.4 to 4.2, the 
star brightened in V by ~ 0.5 m , which is entirely consistent with 
a brightening of a source of continuous emission on the stellar 
surface. 

3.2. Veiling versus line flux and phase 

The equivalent widths (EW) of the emission lines in RU Lup 
vary with time. We have measured EWs for a large number of 
different types of lines formed in different regions around the 
star and investigated how EW depends on VF (as measured in 
the V band). In addition, we have studied how the line inten- 
sities correlate with VF across the lines. The method of cor- 
relation analysis was first described by Gameiro et al. (2006). 
Fig. |2] illustrates these dependancies for two lines, namely He i 
A 5876 A and [S n] A 6731 A. The [S n] line is thought to form 
in an extended wind, where the red-shifted emission is occulted 
by the disk, while the central He i emission is thought to form 
close to shocked regions on the stellar surface. Part of the blue 
wing in He i may trace a central stellar wind. Changes in EW 
are not necessarily related to changes in line flux. In the case 
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Fig. 2. Top panels: Equivalent width (EW) versus veiling fac- 
tor (VF), measured in the V band, for emission lines of He i 
and [S n] for RU Lup, where individual runs are marked. The 
dashed curves correspond to the case of a variable stellar con- 
tinuum but a constant emission line flux. Bottom panels: The 
corresponding correlation functions together with the average 
spectrum (shaded), showing how intensity across the line pro- 
files correlates with veiling. The two lines show opposite corre- 
lations in the blue wings, but only when VF i= 3.5. The narrow 
peaks in the correlation functions reflect the behavior of absorp- 
tion lines not visible in the shaded spectrum. 



of a constant line flux but a variable continuous excess emis- 
sion the EWs would follow the dashed curves in Fig. [2] There 
are long-term changes in the general level of line emission, but 
during states of moderate veiling, VF ^ 3.5, the [S n] lines tend 
to follow the case of constant line flux and varying excess con- 
tinuum. The He i lines show the opposite trend in that during 
states of low to intermediate veiling EW increases with VF. This 
is particularly true for the blue-shifted part of the emission, indi- 
cating that an increase in the accretion rate/accretion luminosity 
leads to increased emission from outflowing He i gas. However, 
when the veiling becomes large EW stays constant with increas- 
ing VF both for He i and [S n]. The Balmer lines and emission 
from metals, like Fe i and Fe rx, also saturate in EW for values 
of VF ^ 3.5. At lower degrees of veiling the metal lines show a 
larger scatter in EW indicating that intrinsic flux changes occur 
on shorter time-scales, which will be discussed elsewhere. 

We conclude that when the veiling increases to very high lev- 
els, the flux in lines formed in very different regions around the 
star do not respond in the same way. This provides an indepen- 
dent illustration of our finding that high veiling is unrelated to 
emission from a shocked region. Both RW Aur A and RU Lup 
show regular radial velocity changes with periods of a few days 
(Paper 1 and 2). We did not find any correlation between de- 
gree of veiling and phase for these stars. Nor did we find any 
correlation between V magnitude and phase from the available 
photometric data banks. 



4. Discussion 

Optical veiling is commonly interpreted as the result of excess 
continuous emission produced in hot areas surrounding the mag- 
netic poles, where accreting gas is violently halted. The photo- 
spheric absorption lines become diluted by this continuum, and 
consequently the degree of dilution, measured by the veiling fac- 
tors, should provide a measure of the strength of the emission 
continuum. We have found that this prediction may hold true in 
our sample of rather "extreme" CTTS, with strong metallic line 
emission, but only for low to moderate values of the veiling fac- 
tors. During periods when the dilution of absorption lines is very 
large, the veiling factors are unrelated to any excess continuous 
emission. This is evidenced from photometric variability versus 
degree of veiling, and also from the behaviour of emission line 
strength versus veiling. 

When the veiling fluctuates around low to moderate val- 
ues, which happens frequently in RU Lup, there can be a pos- 
itive veiling versus brightness correlation in accordance with the 
standard accretion model, and the equivalent widths of emission 
lines vary accordingly. Blue-shifted He i emission becomes en- 
hanced when the excess emission increases, unlike the forbidden 
lines arising in an extended wind. The enhanced He i emission 
indicates that a central stellar wind component is excited during 
an accretion event, while the extended wind component is not 
affected. However, when the veiling increases to large values the 
strengths of all different types of emission lines remain constant. 

In AA Tau, the comparatively weak veiling varies periodi- 
cally with the rotational period (Bouvier et al. 2007b), and the 
line strength of the He i emission correlates with the veiling, sug- 
gesting it originates at a hot surface region moving with rotation. 
In RW Aur and RU Lup, we found dramatic changes in the de- 
gree of veiling even during one night, but no evidence that the 
veiling varies with rotational phase. 

Therefore, "extreme" CTTS may host sources of variable 
continuous excess emission distributed over a small surface as 
predicted in magnetospheric models. However, additional pro- 
cesses are needed to generate the substantial weakening ob- 
served in the photospheric lines. Below we will discuss three 
possible causes. 

A. Variable extinction. If infalling dust grains survive an ac- 
cretion event, then the related increase in extinction in the line- 
of-sight might under specific conditions balance the increase in 
footprint brightness. This possibility was discussed for AA Tau, 
with an almost edge-on disk (Bouvier et al. 120031 |20 07b), and 
for RW Aur A (Paper 2), which on occasion shows strong in- 
verse P Cygni profiles indicative of events of sudden accretion 
in the line-of-sight. However, in RW Aur A there is a lack of de- 
pendence of strength of accretion components in Na D lines on 
veiling (Petrov & Kozack 2007). RU Lup, which is seen nearly 
pole-on (Herczeg et al. 2005), has never shown inverse P Cyg 
profiles. Furthermore, lines formed close to the stellar surface 
and far out in a wind cannot both be occulted by a local accretion 
stream. We find that the case of variable extinction in accretion 
streams is not supported in our study of four CTTS. 

B. Electron scattering. Electron scattering can lead to a di- 
lution of the strengths of photospheric absorption lines without 
contributing excess continuum emission. This mechanism has 
not been invoked before in explaining veiling in T Tauri stars, 
but Stahl & Wolf ( 1980) discussed e-scattering as a cause for 
the broad wings of the Balmer lines in S CrA. The characteristic 
temperature of the infall shock in accreting T Tauri stars is of 
the order of 10 5 K to 10 6 K (see e.g., Calvet & Gullbring [l998l . 
The shocked gas will be highly ionized, so there will be plenty 
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Fig. 3. Spectral differences for different levels of veiling in S 
CrA SE and NW, and in RU Lup. A narrow peak of night sky 
emission enters near 5577 A. In the S CrA stars several photo- 
spheric lines switch from absorption to emission when the veil- 
ing becomes high. Line-filling appears to be present also in cer- 
tain lines in RU Lup, marked with arrows. 



of hot electrons available near the accretion shock. The net ef- 
fect of scattering of light by hot electrons is to shift the affected 
photons to shorter wavelength (the inverse Compton effect) and 
to broaden the frequency distribution (due to the velocity dis- 
persion of the electrons). Code (1949) integrated the radiative 
transfer equation and tabulated the wavelength shift and broad- 
ening of a monochromatic line as a function of electron temper- 
ature. The cross section for electron scattering is the Thompson 
scattering cross section (with minor modifications), so an ion- 
ized mass column of about 2.5 g cirT 2 above the photosphere 
will completely erase all underlying spectral information. Once 
optically thick, electron scattering is a geometric effect, with the 
veiling dependent on the fraction of the surface covered by this 
cloud. As an example, if the continuum excess saturates at VFo 
= 3, one needs to obscure about a fraction l-((l+VFo)/(l+VF)) 
of the star by electron scattering to produce the observed VF. For 
VF = 4, a screen covering 20% of the visible surface is needed. 
Rather special conditions are required for this to occur, such as 
the hot plasma spreading inside the magnetic funnels after an ac- 
cretion event has ceased. In particular, it is difficult to envision 
how the largest values of VF, requiring the largest surface cov- 
erage, can build up. Here, we point out that the process as such 
may be of some importance in TTS, and a more thorough dis- 
cussion is reserved for a forthcoming paper. As described below, 
it appears that yet another process affects the veiling measures, 
and can dominate the dilution of absorption lines in CTTS with 
strong metallic emission lines. 

C. Emission reversals in metallic line cores. In Fig. [3] we 
compare the same spectral region obtained at high and low levels 
of veiling in the S CrA stars and RU Lup. At maximum veiling 
several lines of Fe i are clearly reversed in both S CrA stars, and 
in RU Lup the same lines appear to be filled-in by line emission. 
No clear dependence of line depth on excitation potential was 
found for this star, but strong lines are affected most. We con- 
clude that variable narrow emission components filling-in the 
absorption lines obviously lead to large measured veiling fac- 



tors without any increase in the continuous excess from any hot 
surface area. An increase in line emission will lead to a small 
increase in brightness. In RW Aur A line-filling is not conspicu- 
ous, but may take place (see Fig.Q]). 

The emission line components are narrow, indicating a 
"chromospheric" origin, and that the atmospheric structure 
changes with time. Similar phenomena are seen in other types 
of objects, and in the case of CTTS it was once demonstrated 
by Cram (1979) and Calvet et al. ( 1 19841 ) how unusual spectral 
emission signatures could arise from changes in the T - t re- 
lation in the stellar atmosphere when the temperature minimum 
falls below r = 1 (denoted the "deep chromosphere" by Herbig 
( 11970b ). We found in Sec. 3 that even at high veilings there is 
often a weak correlation between brightness and veiling. These 
variations may relate to minor changes in the accretion rate, 
which in turn could affect the stellar atmosphere surrounding the 
accretion shock leading to the large veiling factors discussed. 
We have found periodic velocity changes in some narrow emis- 
sion components, and the case of variable amounts of chromo- 
spheric emission surrounding the magnetic poles will be further 
explored in a separate paper. 



5. Conclusions 

Our spectroscopic and photometric observations of four clas- 
sical T Tauri stars with very strong emission signatures show 
that the degree of veiling is variable, but when the veiling be- 
comes strong it is only weakly dependent on the stellar bright- 
ness. Measured emission line strengths support this conclusion. 
We demonstrate that veiling can be caused by narrow emission 
lines that fill in the stellar absorption lines and not by a continu- 
ous source of emission, as is usually postulated. We discuss two 
other effects leading to a reduction of absorption line strengths, 
and which may play some role in CTTS. We conclude that the 
degree of veiling, as measured from absorption lines, leads to an 
overestimation of the accretion rates in CTTS with rich emission 
line spectra, especially when the veiling factors VF exceed ~ 3, 
as happens frequently in the stars investigated. In extreme cases 
the overestimation may amount to an order of magnitude. 
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